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Abstract 

Nonequilibrium phase transitions are discussed with emphasis on general features such as the 
role of detailed balance violation in generating effective (long-range) interactions, the importance of 
dynamical anisotropies, the connection between various mechanisms generating power-law correla- 
tions, and the emergence of universal distribution functions for macroscopic quantities. Quantum 
spin chains are also discussed in order to demonstrate how to construct steady-states carrying 
fluxes in quantum systems, and to explain how the fluxes may generate power-law correlations. 
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I. INTRODUCTION 



There are many reasons for studying nonequilibrium phase transitions. Let us start be 
mentioning a few which carry some generahty. 

First and foremost, equihbrium in nature is more of an exception than the rule, and 
structural changes (which constitute a significant portion of interesting phenomena) usually 
take place in nonequilibrium conditions. Thus there is much to be learned about the complex 
ordering phenomena occurring far from equilibrium. 

Second, while very little is understood about the general aspects of nonequilibrium sys- 
tems, the equilibrium critical phenomena have been much studied and have been shown to 
display universal features. This universality emerges from large-scale fluctuations in such 
a robust way that one can expect that similar mechanisms will work in nonequilibrium 
situations as well. Thus, investigating the similarities and differences of equilibrium and 
nonequilibrium orderings may help to discover the distinguishing but still robust properties 
of nonequilibrium systems. 

Third, power law correlations are present in many nonequilibrium phenomena and there 
have been many attempts to explain these correlations through general mechanisms. Closer 
examination, however, usually reveals a close connection to equilibrium or nonequilibrium 
critical phenomena. 

My lectures were designed to revolve around problems related to the above points. The 
lectures are built on the theory of equilibrium phase transitions 0], thus I assume knowl- 
edge about both static and dynamic critical phenomena at least on the level of familiarity 
with the basic concepts (symmetry breaking, order parameter, diverging correlation length, 
order parameter, scaling, universality classes, critical slowing down, dynamical symmetries). 
I discuss simple examples throughout so that enterprising students could try out their luck 
implementing their ideas in simple calculations. Due to space restrictions, however, not all 
of the details discussed in the lectures and afterwards are covered in the written notes. In 
particular, the picturesque parts of the explanations are often left out since they take up 
disproportionally large part of the allowed space. Nevertheless, these pictures may be im- 
portant in both understanding and memorizing, and I strongly encourage the reader to go 
through the slides of the lectures, as well. They can be found through the homepage of the 
school, |http://dpm.univ-lyonl.fr/houches_ete/lectures/| or at |http:/ /poe. elte.hu/~racz/ . 



A. Nonequilibrium steady states 

A general feature that distinguishes a nonequihbrium steady state (NESS) from an equi- 
hbrium one is the presence of fluxes of physical quantities such as energy, mass, charge, etc. 
Thus the study of NESS is, in a sense, a study of the effects of fluxes imposed on the system 
either by boundary conditions, or by bulk driving fields, or by some combination of them. 
A nonequilibrium steady stateswell known example is shown on Fig.|I|. 

heat bath T2 



bulk field 



viscous fluid 




NESS 



heat bath Tj 
FIG. 1: Setup for Rayleigh-Benard experiments. 

This is the Rayleigh-Benard experiment in which a horizontal layer of viscous fluid (the 
system) is heated from below i.e. it has two heat baths of temperatures Ti and T2 attached 
(boundary conditions generating an energy flux). The presence of gravity (the bulk drive) 
is also important (it generates mass and momentum fluxes at large 6T = Ti — T2 > 0). 

For Ti = T2 this system relaxes into a quiescent equilibrium state while a small ST will 
also result in a quiescent state but it is already a NESS since energy flux is flowing through 
the system. Increasing 6T, this steady state displays a nonequilibrium phase transition 
(Rayleigh-Benard instability), first to a stationary convective pattern, and then to a series 
of more complicated structures which have fascinated researchers for the past century 

Starting from the Navier-Stokes equations, one can arrive at a mean-field level of under- 
standing of the above phenomena. It is, however, not the level of sophistication one got 
used to in connection with equilibrium phase transitions. There, we have simple exactly 
solvable models such as e.g. the Ising model which give much insight into the mechanism 
of ordering and, furthermore, this insight can be used to develop theories which reveal the 
universal features of equilibrium orderings 0. 

The trouble with the Rayleigh-Benard system is that we do not have a theory even 
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for the NESS. The reason for this is that the fluxes result in steady state distributions 
, P*, which break the detailed balance condition Wn-^n'Pn = Wn'-^nPn'^ where n and n' 
are two "microstates" and Wn-,n' is the rate of the n ^ n' transition. As a consequence 
of the breaking of detailed balance, a NESS is characterized not only by the probability 
distribution, P*, but also by the probability currents in the phase space. Unfortunately, we 
have not learned yet how to handle the presence of such loops of probability currents. 

The main lesson we should learn from the Rayleigh-Benard example is that, in order to 
have Ising type models for describing phase transitions in NESS, one should use models 
which relax to steady states with fluxes present. Such models have been developed during 
the last 20 years, and most of my lectures are about these stochastic models defined through 
"microscopic" elementary processes. The first level of description is in terms of master equa- 
tions which are conceptually simple and allow one to make use of general results (uniqueness 
of stationary state, etc.) which in turn are helpful in defining dynamics that leads to NESS 
(Sec.0). The next level is to describe the same problem in terms of Langevin equations 
and develop field-theoretic techniques for the solution. Our discussions will include both 
levels of description and I hope that at the end an understanding will emerge about a few 
results which grew in importance in the last decade (generation of long-range interactions 
and effects of dynamical anisotropics (Sec.y), connection between mechanisms generating 
power-law correlations (Sec.^), and universality of distribution functions for macroscopic 
quantities (Sec.[rVD). 

Before starting, however, I would like to insert here a little essay about effective temper- 
atures. This concept is being widely discussed in connection with slowly relaxing systems, 
the topic of this school. So it may be of interest to present here a view from the perspective 
of NESS. 

B. Problems with usual thermodynamic concepts 

Systems close to equilibrium may retain many properties of an equilibrium state with 
the slight complication that the intensive thermodynamic variables (temperature, chemical 
potential, etc.) become inhomogeneous on long lengthscales and they may slowly vary in 
time. This type of situations are successfully dealt with using the so called local equilib- 
rium approximation ||^, with the name giving away the essence of the approximation. The 
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applicability of the concept of local equilibrium should diminish, however, as a system is 
driven far from equilibrium. Nevertheless, questions of "how large drive produces a far- 
enough state" and "couldn't one try to find a new equilibrium state near-by" are regularly 
asked and have legitimacy. So I will try to illuminate the problems on the example of the 
fluctuation-dissipation theorem much discussed nowadays due to attempts of associating ef- 
fective temperatures with the various stages of relaxation in glasses or with steady states 
in granular materials 

Let us consider a simple system of Ising variables a with Hamiltonian Tig and in equi- 
librium at temperature f3 = 1/(A;bT). Assuming that there is an external field H coupling 
linearly to the macroscopic magnetization M = cTj, one can write the equilibrium distri- 
bution function as 

Pe,ia) = Z-h-P^'^^^^+^^^'^^^^ (1) 
The average value of the magnetization is given by 

(M) = M(a)e-'^^"('^)+^^'^('^) . (2) 

and the static limit of the fluctuation-dissipation theorem is obtained as 

d{M) 



OH 



= m') (3) 



where we assumed the system to be in the high-temperature phase ((M) = 0). Note the 
simplicity and the accompanying generality of this derivation. It uses only the fact that the 
external field is linearly coupled to the quantity (M) we are considering. 

The fluctuation-dissipation theorem is used in many ways. It helps simplify field-theoretic 
studies of fluctuations through diagrammatic expansions and it also gives a powerful checking 
procedure in both experiments and Monte Carlo simulations. Note that eq.(^ can also be 
used to define the temperature of the system through (3 = Xm/{^'^)i ^"^^ temperature 
defined in this way would be the same when using different "M"-s and conjugate fields "if". 

It is clear that a fluctuation-dissipation theorem generalized to NESS would be extremely 
useful. Let us now try to imagine how a similar relationship may arise when we drive 
the above system away from equilibrium (e.g. by attaching two heat baths of different 
temperatures). If the system relaxes to a NESS then there will be steady-state distribution 
function, (P) but the effective Hamiltonian (InP) will contain all the interactions allowed 
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by the symmetries of the system. Thus, assuming that the effective Hamiltonian can be 
expanded in if, one finds in the H —>■ hmit 

where 5*3 is a notation for sums over all three-spin clusters with different couplings for 
different types of spatial arrangements of the three spins. Furthermore, a, b and all other 
newly generated couplings depend on the original couplings in Ho and on the temperatures 
of the heat baths. 

Now a derivation of the fluctuation-dissipation theorem similar to the equilibrium case 
yields a more complicated equation 

d{M) 



dH 



= a[(M2) + (M^3) + ...] . (5) 



There are two ways a simple form for the fluctuation-dissipation theorem may emerge from 
eq.(^). One is that a nonlinear field 

Q = M + S^ + ... (6) 
that is conjugate to H can be introduced (and effectively worked with). Then one obtains 

Xa = «(^?') (7) 

and thus a becomes the nonequilibrium (3. 

The other possibility is that a mean-field type decoupling scheme works well and then 

{MS,) = {M')fiC2) (8) 

where /(C2) is a functional of the two-point correlations (and similar expressions are ob- 
tained for {MS2n+i))- Then equation (|^) becomes 

X,, = aJ^{C2){M'). (9) 

If the theory provides T{C2) then the effective temperature can again be read of from the 
above generalized fluctuation-dissipation relationship. 

There are problems with both lines of reasoning. Apart from the practical difficulties 
of nonlinear fields and the validity of mean-field type approaches, the main conceptual 
difficulty is the fact that changing from the magnetization to other fields (e.g. energy) the 



nonequilibrium version of the fluctuation-dissipation theorem leads to different values for 
the same "temperature" ||T0|. There are cases where the above schemes generalized to time- 
dependent processes work both at the theoretical 0, H and the experimental levels 
but there are clear examples when the concept of effective temperature does not apply []TU 



Thus the meaning and use of nonequilibrium temperature has not been clarified enough to 
make a verdict on it. 



II. PHASE TRANSITIONS FAR FROM EQUILIBRIUM 

As mentioned in the Introduction, nonequilibrium phase changes constitute a large part 
of interesting natural phenomena and they are studied without worries about wider con- 
texts. From a general perspective, on the other hand, the investigations of nonequilibrium 



phase transitions |Tl|, can be viewed as an attempt to understand the robust features 
of NESS. This view is based on the expectation that the universality displayed in critical 
phase transitions carries over to criticality in NESS as well. If this is true then studies of 
the similarities to and differences from equilibrium will lead to a better understanding of 
the role and general consequences of the dynamics generating NESS. 

In the following subsections, we shall construct, describe, and discuss models which dis- 
play nonequilibrium phase transitions. Apart from getting familiar with a few interesting 
phenomena, the main general conclusion of these discussions should be that dynamical 



anisotropics often yield dipole-like effective interactions [|13t 0, |1^ and, furthermore, com- 
peting non-local dynamics (anomalous diffusion) generates long-range, power-law effective 
interactions [|16|. Along the way, we shall also understand that the detailed-balance violating 



aspects of local relaxational dynamics do not affect the universality class of the nonequilib- 
rium phase transitions |]15|, p!7|] . 



A. Differences from equilibrium - constructing models with NESS 

The violation of detailed balance has the consequence that not only the interactions 
determine the properties of the NESS but the dynamics also plays an important role. In 
order to understand and characterize the role of dynamics, a series of simple examples will 
be discussed in the following subsections. 
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First, let us discuss how to construct a model which yields a NESS in the long-time limit. 
A simple way is to attach two heat baths to a system, each generating a detailed-balance 
dynamics but at different temperatures. To see an actual implementation, let us consider 
how this is done for the one dimensional kinetic Ising model. This type of models have been 
much studied and a collection of mini-reviews about them can be found in [113 . 



The state of the system {cr} = {..., cxj, o"j+i, . . .} is specified by stochastic Ising variables 
crj(t) = ±1 assigned to lattice sites i = 1,2, . . . , N. The interaction is short ranged (nearest 
neighbor) —Jdiai^i and periodic boundary conditions ((Tat+i = ai) are usually assumed. The 
dynamics of the system is generated by two heat baths (labeled by a = 1, 2) at temperatures 
Tq, meaning that the heat baths try to bring the system to equilibrium at temperature Tq, 
by e.g. spin flips and spin exchanges, respectively. 

Let us denote the rate of the flip of i-th spin {at — > —<Ji) by wl^\{cr}), and let the rate of 
the exchanges of spins at sites i and j {ai ^ aj) be w['^\{a}). Then the dynamics is defined 
by the following master equation for the probability distribution P{{a},t) : 



dtP{W},t) =E. wl'\{ah) P{{ah,t) -w','>{{a}) P{{a},t) 



(1) 



+ 



,(2) 



(2), 



E., K'iWh) PiWh,t) - <({a}) P{{a},t) 



(10) 



where the states {cr}i and differ from {cr} by the flip of the i-th spin and by the 

exchange of the i-th and j-th spins, respectively. 

The assumption that the dynamics is generated by heath baths means that the rates 
satisfy detailed balance at the appropriate temperatures: 



(11) 



where P^'^ ~ exp[—J/TaJ2i^i^i+i] the equilibrium distribution of the Ising model at 
temperature Tq,. Eq.(|TT]) leaves some freedom in the choice of w^-s, and one is usually 
guided by simplicity. The most general spin flip rate that depends only on neighboring spins 



has the following form |T9 



1 
2^ 



(12) 



Without any other heath baths, equations (|T0D and (0) deflne the Glauber model |T9 



which relaxes to the equilibrium state of the Ising model at temperature Ti deflned through 
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7 = tanh(2J/A;BTi). The time-scale for flips is set by ti and S is restricted to the interval 
-1< 5 < 1. 

The competing dynamical process is the generation of spin exchanges (Kawasaki dynamics 
pip by a second heath bath at a temperature T2 7^ Ti. In the simplest case, the exchanges 
are between nearest neighbor sites and the rate of exchange satisfying detailed balance (|TTp 
is given by 

V) = ^ ^ ~ y + (^i+i(^z+2) ■ (13) 

where 72 = tanh(2J/fcBT2) and T2 sets the timescale of the exchanges. It is often assumed 
that the exchanges are random (T2 = 00) and thus w^'^\a) = l/(2r2). 

Equations (|IOD , (|1^) and ([T3|) deflne a model which can be shown to have a NESS and 
one can start to ask questions about the phase transitions in this steady state. The gen- 
eralizations to higher dimensions, to various combinations competing dynamics (flip - flip, 
flip - exchange, exchange - exchange), and other types of dynamical steps (resulting e.g. 
from a bulk driving fleld) should be obvious. Just as it should be obvious that there are 
not too many exactly solvable models in this fleld and most of the results are coming from 
simulations [|ll], [T^ . 

Before turning to results, let us also introduce a Langevin equation description of compet- 
ing dynamics. The Langevin approach has been successful in dynamic critical phenomena 
where the counterparts of the Glauber and Kawasaki models are called Model A and B, 
correspondingly. This correspondence makes the "two-heath-baths" generalization straight- 
forward. The coarse grained magnetization of the Ising model is replaced by ra-component 
order-parameter fleld, S'*(x, t), {i = 1, ...,n) (n = 1 is the Ising model, and n — 00 is the 
spherical limit that allows simple analytic calculations as shown below). The system evolves 
under the combined action of local relaxation (Model A) satisfying detailed balance at tem- 
perature Ti and diffusive dynamics (Model B) at temperature T2 which yields the following 
Langevin equation for the Fourier transform Sl^{t) 

9i^; = -£«5; + r/i(q,t)-42)^^ + ^^(q,t) . (14) 

Here C'^^ Sl^ = vi^^ SF^"^ /6S^_q^ with F^") being the free energy at temperature Tq,, and 
the kinetic coefficient Fq"'' is enforcing the conservation laws. In particular, Fq"^ = Fq°^ in 
Model A without conservation laws, and Fq"'' = D^'^^c^^ for Model B with diffusive dynamics 
conserving the total magnetization. In case F*^"^ is the Landau-Ginsburg functional, we have 
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q' q" 



(15) 



where Tq is linear in Tq and m ~ l/n in the n —>■ oo spherical limit. In order to ensure that 
in case of a single heat bath, the system relaxes to equilibrium satisfying detailed balance, 



the noise terms in eq.(|T4| are Gaussian-Markovian random forces with correlations of the 
form: 

R(q, tWA^, t')) = 2r(f )<5„„,5.,<5(q + q')5(t - t') . (16) 

Eqs. (p!4D,(p!5D and (|T6|) define the model for the particular competing dynamics chosen 
and we are now ready to deduce some features of the NESS generated. Of course, just as in 
case of kinetic Ising models, the number of possible competing dynamics is infinite and the 
question is whether conclusions of some generahty could be reached. 



B. Generation of long-range interactions - nonlocal dynamics 

The remarkable consequences of competing dynamics can be seen already on the example 
of d = 1 flip-and-exchange model which may produce ordering even though the interactions 
are of short range. Indeed, if Ti temperature spin flips are competing with T2 = oo spin 
exchanges of randomly chosen pairs then the system orders below a certain Tic 



It turns out that the transition is of mean-field type and this gives a clue to understanding. 
Indeed, let us imagine that the rate of spin exchanges is large compared to the rate of flips. 
Then the random exchanges mix the spins in between two flips and the flipping spin sees 
the "average spins" in its neighborhood - a condition for mean-fleld to apply. 

The mean fleld result can also be interpreted as the generation of inflnite-range effective 
interactions. This interpretation can be put on more solid base by studying the above model 
with T2 = 00 spin exchanges where the probability of exchange at a distance r is decaying 
with r as p{r) ~ i^^d+a ^^j^^ spins exchanges are a dimensional Levy flights in dimension d). 
The system orders again below a Ti^ and the examination of the critical exponents reveals 
|T6| that the transition is in the universality class of long-range interactions decaying with r 



as J(r) ~ l/r*^"*""". It is important to note that the above results are nonequilibrium effects 
which would disappear if the spin exchanges would also be at Ti. 
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Let us now see if the same results can be derived from the Langevin equation approach. 
The spin flips are translated into the Model A part of the dynamics while the Levy flights 
can be represented by anomalous diffusion with Fq •* = D^'^^c^^ replaced by Fq = D^'^^q^ 



with < cr < 2 being the dimension of the Levy flight. Thus the Langevin equation becomes 



— -^OV'U^yy^q ^0"^ J J ^q'^q"^q-q'-q' 

q' q" 

-Dq-Si + f^{t) . (17) 

Note that due to the randomness of the Levy flights (T2 = 00), the interaction and the 
nonlinear terms are missing in the Levy flight part (second line) of the equation. As discussed 
in Sec. p A| , the rj-s are Gaussian-Markoffian random forces with correlations ('7q(t)^7q/(i')) — 



2Fo5.,5(q + q')5(t - t') and «(t)r/^,(t')) = 2Dq''6,,6{q + q')5(t - t'). 

In order to see the generation of long range interaction in the above model, let us first make 
an exact calculation |TB| in the spherical limit (n 00) where fluctuations in m ^ S^(t)S^^,(t) 



may be neglected and this quantity may be replaced by 

n 

^J2{Si^{t)S'^,{t))=unC{q,t)6{ci + c^) , (18) 



n 

U 



where the brackets () denote averaging over both the initial conditions and the noises 77 
and f] (we restrict ourselves to the study of the high-temperature phase where the dynamic 
structure factor C{q,t) = {S^it)S-'_q^(t)) is independent of j). 

The decoupling (|T^) leads to a linear equation of motion and so the self-consistency 
equation for C (g, t) can be easily derived 

C{q, t) = 2(Fo + Dq"") / dt'e'^^'^ {To[roW+unC{q,s)]+Dq'^]ds _ ^^9^ 

Jo 

Here, the initial condition C(g, 0) = was used for simplicity. The t ^ 00 limit does 
not depend on the initial condition and the equation for the steady state structure factor 
C{q) = C{q,t ^ 00) becomes 

r( ] = To + Dg- 

Fo(ro + q^ + unS) + Dq- ' ^ ' 

where S = f dqC{q). 

The long- wavelength instabilities are determined by the g form of C{q) which for 
< 0" < 2 can be written as 

C{q)^{ro + Xq'' + unS)-^ , (7) 
12 



with A = D/Tq. This form coincides with the long- wavelength limit of the equilibrium 
structure factor of a spherical model in which the interactions decay with distance as r~'^~'^. 
Consequently, both the self-consistency equation for r = tq + unS and the critical behavior 
that follows from it are identical to that of the equilibrium long-range model. Thus we can 
conclude that the critical properties of the NESS are dominated by an effective long-range 
potential proportional to r~'^~°". 

The above conclusion should be valid quite generally for finite n as well. Looking at 
6q-(0)' ^^"^ ihaX the correlations in the effective noise {f]eff = f] + fj) have an 
amplitude 2(ro + Dq'^). One expects that the Dq'^ term can be neglected in the long- 
wavelength limit and thus that the noise fj in the Levy-flight exchanges can be omitted. 
Without fj, however, the system described by equation (2) satisfies detailed balance and has 
an effective Hamiltonian which, apart from the usual short-range interaction pieces, contains 
the expected long-range part A / c?qg'^^^5'*(q)S''(— q). 

We should note here that a — 2 corresponds to usual diffusion and that the above 
arguments changes for cr = 2 since no long-range interactions are generated any more, and 
no change in critical behavior occurs. This result is another way of saying that Model A type 
dynamics is robust against diffusive perturbations which break detailed balance Note 
also that if both the competing dynamics are relaxational then, adding up the corresponding 
deterministic and noise parts in the Langevin equation (|Tip, one can easily deduce that the 
breaking of the detailed balance does not change the universality class of the equilibrium 
phase transition. 



C. Generation of long-range interactions - dynamical anisotropies 

In order to understand the meaning of dynamical anisotropy, let us consider the two- 
temperature, diffusive kinetic Ising model on a square lattice. Two heat baths are 
attached and both of them generates nearest neighbor spin exchanges. Exchanges along 
one of the axis (called 'parallel' direction) satisfy detailed balance at temperature Ty while 
exchanges in the 'perpendicular' direction are produced by a heat bath of temperature T^. 
It is important to note that the interactions JaiCj are the same along both axes. It is the 
dynamics that is anisotropic. 



For T± = Til = T, this is the Kawasaki model |20| which relaxes to the equilibrium 
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Ising model at T and, consequently, it displays a continuous transition. Since the dynamics 
conserves the total magnetization, the ordering for T < Tc appears as a phase separation. 

For Tj_ 7^ T||, on the other hand, there is a flow of energy between the || and _L heath 
baths and the system relaxes to a NESS. MC simulations |T^, ^ show that a critical phase 



transition is present for Tj_ 7^ Ty as well but the phase separation is distinct from that 
occurring in equilibrium. The interfaces between the domains of up and down spins align 
with normals along the directions of lower temperatures. Thus the symmetries of the ordered 
states are different from the symmetry of the equilibrium order where interfaces with normals 
along any of axes coexist (isotropic ordering). As a consequence, the universality classes of 
the II and _L orderings are found to be distinct from the Ising class Renormalization 
group calculations actually show that the universality class of the nonequilibrium ordering 



is that of a uniaxial ferromagnet with dipolar interactions |24|. 

A dramatic demonstration of the long-range nature of the interactions generated by 
anisotropic dynamics comes from the generalization of the above model to n = 2 component 
spins (2T — XY model). One finds that the NESS in this system displays an ordering 
transition [1^, a fact that would be in contradiction with the Mermin- Wagner theorem |25 
should the effective interaction be short-ranged. 

Let us try now understand the generation of the dipole-like interactions using the 
Langevin equation approach and considering the spherical limit again. The two-temperature, 
diffusive Ising model corresponds to competition of two Model B type dynamics along || and 
± axes. Thus the equation of motion is 5'q(t) = £||(g)S'q + ?7*(q, t) + £_L(g)S'q + ?7^(q, t) with 
the diffusion in the a = ||, ± directions described by the corresponding Ca terms: 



(r- + q')Sl + ^T.J J 
i = l q/ q// 

2 



(21) 



where q = (gy, gx). Note the isotropy of the interaction (g term) and the different temper- 
atures {tq ) for diffusion in different directions. The noise correlations follow from detailed 
balance requirements, (^^^(q, t)^]^/ (q', t')) 
= 2D^ql5^^.5,j5{ci + q')5(t - t'). 

Just as in the kinetic Ising model for Tq = r^, we have Model B with anisotropic diffusion 
(not a dynamical anisotropy!) and with an equilibrium steady state, while a NESS is pro- 
duced for Tq Tq. The nature of the phase transitions in the NESS becomes transparent in 
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the spherical limit (n — oo and m ~ 1/n) where the fluctuations in w ^ 5'^(t)S'q,(t) may be 
neglected. This linearizes the equation of motion and allows to write down a selfconsistency 
equation for C(q) = {S^{t)S^^,(t)). The t —>■ oo limit then yields the steady state structure 
factor in the following form [|I^ 

C(q) = — » ^ , (22) 

where a = D±/D^^ and S = un J (iqC(q). 

One can see now the origin of dipole-like effective interactions. Because of the dynamical 
anisotropy, the q ^ limit of C(q) is different whether flrst q\\ —>■ and then q±_ — *■ 
or vice versa. This singularity of the long-wavelength limit translates into such power 
law correlations in real space which are characteristic of dipole interactions in equilibrium 
systems. Hence the conclusion [jl^, ^ that the dynamical anisotropy has generated dipole- 
like interactions. Note that this is a nonequilibrium effect. The long- wavelength singularity 
disappears as soon as the heat baths have equal temperatures {r^ = r^). 

The dynamical anisotropy is a strong effect and its mechanism of action is rather simple 
as we have seen above. Accordingly, it is the most viable candidate to change the universality 



class of equilibrium phase transitions by breaking detailed balance |15 



D. Driven lattice gases, surface growth 

Models of NESS have a long history but the flrst one that became the center of attention 



and was recognized as the "Ising model" of NESS was the driven lattice gas (see |TT[ and 
references therein). The model can be understood as a kinetic Ising model with the up-spins 
being the particles in the lattice gas. Spin exchange dynamics at temperature T represents 
the particle diffusion and an external bulk fleld Er^ drives the up-spins (particles) along one 
of the lattice ). In order to have a NESS with particle current one must also use 

periodic boundary conditions in the fleld direction. 

This model displays a critical phase transition in its NESS and the phase separation that 
follows at low T is characterized by strong anisotropy: the interfaces align parallel with E^. 



Thus one can see some similarities with the two-temperature diffusive model of Sec. [II C 



One can easily recognize that dynamical anisotropy is at work here. The driving fleld can 
be considered as a second heath bath which generates the essential part of the dynamics along 
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the X direction. There is a difference, however, from the two-temperature model of Sec. [II C 
in that the drive now has a directionahty (the forward-backward symmetry of diffusion is 
broken). As a resuh the phase transition is expected to belong to a new (nonequilibrium) 
universality class distinct from that of the dipole class. This is indeed what has been 
obtained in renormalization group calculations [|n|]. Unfortunately, the structure of the 
long- wavelength singularities in the two systems are similar and thus there are difficulties in 
observing the differences in numerical work. This problem has generated some debate that 



is still going on |2g]. We believe the debate will not modify the general picture summarized 
in 



Tl[| , and it will not change the conclusion about the importance of dynamical anisotropy. 



An interesting and important field where even the simplest systems show "effective" 
critical behavior due to the unbounded long-wavelenth fluctuations is the field of surface 
growth processes. Most of the roughening transitions and transitions between various rough 



phases are genuine nonequilibrium phase transitions and have been much studied [27, 28 



Remarkably, however, these transitions have not provided new insight into the general fea- 
tures of nonequilibrium criticality, they merely confirmed that the dynamics and dynamical 
anisotropy play an important role in determining the universality classes of growth processes. 

One of the unsolved problems that continues to fascinate researchers is the phase tran- 
sition in the Kardar-Parisi- Zhang equation We shall discuss the problem of growing 
surfaces including KPZ equation in connection with the nonequilibrium distribution func- 
tions in Sec.|IV|. 



E. Flocking behavior 

Up to this point, we have considered usual physical systems driven out of equilibrium. 
Here I would like to give a taste of what awaits one if the studies are extended to the living 
realm. 

Living creatures can also be viewed as units attached to two heat baths. One of them is 
the internal energy source which on a short time-scale is an infinite bath from which energy 
can be drawn at a given rate. The other one is the surroundings to which energy is lost by 
dissipation (friction, heat loss, ...) due to the activity of the unit. This view suggests that a 
collection of such self-propelled units will show orderings (nonequilibrium phase transitions) 
depending on the interactions between the units, on their density, on their possible motions 
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and on the dissipation mechanisms. Indeed, collective behavior is often observed in flocks 
of birds, in schools of flsh, in swimming cells, etc. and, as shown below, some of these 
phenomena can be described in terms of a surprisingly simple model. 

Model PD| was introduced to describe the collective motion of self-propelled particles 
with birds and bacteria being candidates for these particles. We shall use the language of 
"birds" below. 




FIG. 2: Flocking: Trajectories of 10000 self-propelled particles in the model described in the text. 
The parameters are chosen so that the stationary order parameter is 4> = 0.8. Each particle is 
represented by a point marking its current position as well as a continuous line showing its recent 
(10 time step long) trajectory, (a) Initial stage of the relaxation, (b) the stationary regime. Pictures 
are curtesy of A. Czirok and T. Vicsek. 

The basic assumptions of the model are that (i) the birds fly with constant speed liTjl = 1 
{i is the bird index), and (ii) the birds adjust their direction 6i in time intervals of r = 1 to 
the average direction of other birds within a distance r 

e,(t + T) = {e{t))r + vi (23) 

where rji is random noise with amplitude rj. 

Assumption (i) handles the energy in- and outflow by strictly equating them, while as- 
sumption (ii) handles the interactions by seemingly reducing them to interactions in the 
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space of velocity directions. This is not quite so, however, since the motion of the birds 
Xi{t + r) = Xi{t) + Vi{t)T couples the directional and spatial motions. 

The control parameters in the system are r, 77 and the density of particles p. Keeping r 
and p fixed while varying the "temperature" t], one finds that the birds are fiying randomly 
for 7] > ric{r, p) while collective motion develops below t]c where the birds tend to move in the 
same direction. An order parameter characterizing this spontaneously symmetry breaking 
can be chosen e.g. (f) = \ Vi\/N. Fig.^ shows the time evolution deep in the ordered 
regime (0 = 0.8) starting from a random configuration. One can see that local orientational 
order develops in the initial stages of relaxation (the state here shows resemblance to the 
states in classical XY ferromagnet) while the stationary state with almost full orientational 
order shows large density fiuctuations. The structure and the large density fiuctuations 
observed in the ordered state and, furthermore, the measurements of the critical exponents 



of the transition ||3^ suggest that the ordering in this system is in a universality so far not 
encountered. 

A remarkable field theory has also been constructed for fiocking It is a generalized 
Navier-Stokes equation with additional Model A type terms which drive the velocity to 
{\vi\) = 1. This theory explains the large density fiuctuations present in the ordered state. 
The investigation of ordering transition is at a higher level of difficulty, however, and has 
not been completed yet. 

Clearly, much remains to be done before we understand fiocking and before the model 
can be compared with experiments quantitatively. Nevertheless, activity is expected in 
this direction since the model of fiocking is not much more complicated than the more 
standard NESS models discussed above, and, at the same time, it has close connection with 
experimentally observable, truly " far-from-equilibrium" phenomena. Hopefully, by designing 
and understanding similar models, a kind of " universality map" of the collective dynamics 
of self-driven units can be found. 



III. WHERE DO THE POWER-LAWS COME FROM? 

Systems displaying power law behavior in their various characteristics (correlation in 
space or time, fiuctuation power spectra, size-distributions, etc.) are abundant in nature. 
The most impressive examples are found in biology (e.g. the metabolic rate vs. mass rela- 
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tionship for living creatures displays scaling over 28 decades p^) but there are remarkable 
examples in solid state physics (power spectra of voltage fluctuations when a current is flow- 
ing through a resistor - 6 decades of scaling), in geology (the number of earthquakes vs. 
their magnitude |^ - 5 decades) and scaling over 2-3 decades is seen everywhere [see e.g. 



the white-dwarf light emission the flow of sand through an hourglass the number 
of daily trades in the stock market , water flow fluctuations of rivers PB[ , the spike trains 



of nerve cells |^9|, the traffic flow on a highway interface fluctuations dissipation 



in the turbulent systems 

Understanding the (possibly) common origins of scaling in the above phenomena appears 
to be a highly nontrivial task. Power laws, of course, arise naturally in critical phenomena 
and we understand them: their origins are in the diverging fluctuations at the critical point. 
Thus the first question one may ask is the following. 

• Can the power laws just be the result of nonequilibrium phase transitions and the 
associated critical behavior? 

In equilibrium systems, however, one must tune a parameter to its critical value in order to 
observe scale-invariant behavior while nonequilibrium systems appear to be in scale- invariant 
states without any tuning. Thus the answer to the first question appears to be negative. 

The wide variety of the phenomena in the above list suggests that the next question could 
be as follows. 

• Can the scaling merely be a natural outcome of complex dynamics? 

After all, we have seen in Sec. [1 C that competing dynamics may generate long-range (power- 



law) interactions which may be at the origin of scaling even away from a critical point. The 
answer to this question may be a yes but, unfortunately, many of the problems mentioned are 
not amenable to an analysis in terms of simple competing dynamics and then the following 
question remains unanswered: 

• What are the ingredients of complex dynamics which determine the existence and the 
characteristics (e.g. exponents) of the power laws? 

There is an attempt to answer all the above question in affirmative along a logic that 
begins with the notion of self-organized criticality (SOC) introduced by Bak, Tang and 
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Wiesenfeld |^ . According to this notion, systems with complex dynamics tune themselves 
to a state with a kind of avalanche type dynamics that is underlying a large number of 
scale-invariant phenomena. The notion of SOC has now been understood in terms of an 
interplay between local and non-local dynamics pSf which indeed tunes the system MM to 
a nonequilibrium (absorbing-state) critical point. Then the problem of SOC is reduced to 
investigating the absorbing state transitions P3| and the characteristics of power laws can be 
determined by studying the universality of absorbing state transition. This is an interesting 
and active field of research and it is worth understanding the main points. Accordingly, I 
will discuss SOC in the next subsection, and will explain the connection to absorbing state 
transitions in the following one. 



A. Self- organized criticality (SOC) 

The first model of SOC was introduced to describe sandpiles. Later developments, how- 
ever, made energy packets from the grains of sand, so the balls in Fig.|^ will be grains at 
the beginning but will be called energy packets later. The dynamics defining the model 
consists of local and non-local elements. The sites of a (usually) two-dimensional lattice 
are occupied by grains and the local aspect of their dynamics is in the redistribution of the 
grains. If a site contains more than Zc grains (e.g. Zc = 4 on a square lattice) then the site 
is active and Zc grains are redistributed to the neighboring lattice sites. The redistribution 
of particles (avalanche) continues until active sites are found. Clearly, an avalanche stops 
after a while since the redistribution leads to loss of particles at the boundaries (or, in terms 
of the energy model, dissipation occurs at the boundaries). Once the avalanche stopped, an 
external supervisor notices it (this is the nonlocal part of the dynamics) and starts to add 
new particles (energy is injected into the system) until a new avalanche starts. 

The above dynamics yields a stationary state in the long-time limit, and the steady state 
characteristics of avalanches can be measured. Such a characteristics is e.g. the number of 
sites s which become active during the process, and the remarkable feature of this model is 
that the distribution of s (and of other quantities such as the spatial size and the lifetime 
of the avalanches) is found to display a power law form 

P{s) ~ . (24) 
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FIG. 3: Sandpile model. Particles (energy packets) are deposited on a two-dimensional substrate 
(one-dimensional section is shown). Injection stops when a site becomes active, i.e. it is occupied by 
more than Zc particles. Then redistribution to neighboring sites take place and particles disappear 
(dissipation of energy) at the edges. The process continues until all active sites are eliminated. 
Then the particles source is switched on again. 

Thus one discovers that although the model does not contain parameters to tune, nev- 
ertheless it shows critical behavior {zc can be changed without changing the critical! ty of 
the outcome). This observation generated a large amount of activity and criticality-without- 
tuning was seen in a number of similar models The resulting notion of self- organized 



criticality grew in importance ^^d, accordingly, new effort was put into understanding 
how SOC works. 



An important feature that was recognized quite early is the existence of a non-local 
supervisor who watches the activity of the avalanches and, upon ceasing of the activity, 
switches on the source of particles (or of energy). In principle, non-local dynamics can 
generate long-range correlations in both time and space so the emergence of criticality is 
not necessary a surprise. Viewing the problem from another angle, the non-local dynamics 
separates the timescales of the avalanches and of the particle injection. Thus, in practice, 
there is tuning. Namely, the system is considered in the limit of particle injection rate going 
to zero (actually, the dissipation is also tuned to zero since the particles disappear only at 
the boundaries of the system). 

The zero injection rate, however, does not have to be a critical point, and the next impor- 



tant development was the demonstration that it is indeed a nonequilibrium (absorbing 
state) critical point. 
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B. Absorbing state transitions and their connection to SOC 



Absorbing state transitions appear in many contexts in nonequilibrium statistical physics 
ir^ , |8|, and they are studied intensively since they are thought to be one of the truly 
nonequilibrium phenomena without counterpart in equilibrium systems. In order to under- 



stand the basics of it, let us consider a fixed energy sandpile model |E9[ shown on Fig.|. 



E<Ej, dying out activity 
absorbing state 



E>E^ active phase 



• • • • 



o 



FIG. 4: A one-dimensional fixed energy sandpile model. Dynamics is defined by the energy (par- 
ticles) being redistributed if a site contains more than units of energy (the units which are 
redistributed in the next time-step are marked by large dots in their center). The total energy of 
the system is conserved since the boundaries prevent the loss of energy. 

This model differs from the sandpile model by calling the particles energy units, and by the 
absence of both the injection and the dissipation of energy (particles). Thus the total energy 
E is conserved and as one can easily see from Fig.|^ the behavior of the system is essentially 
different at small and large values of E. At small E < E^, the activity (redistribution) 
ceases in the long-time limit and the system falls into a so called absorbing state. For large 
E > Ec, on the other hand, there are always active sites and the redistribution continues 
forever. For t ^ oo, the system settles into a steady state which is called active state and 
the activity can be quantified by measuring e.g. the number of active sites. One finds then 
that the absorbing state transition (i.e. the absorbing-active state transition) is a critical 
phase transition with the activity changing continously through the transition point E^ (see 
Fig.|). 
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FIG. 5: Activity as a function of the total energy for the fixed energy sandpile model described 
on Fig.^. The evolution of the system as dissipation at the boundaries or the energy injection is 
switched given by the left and right arrows, respectively. 

Once the absorbing state transition is understood, it is easy to make the connection 
to SOC. Indeed, let us assume that we have the fixed energy sandpile in the active state 
{E > Ec) and let us switch on the dissipation at the boundaries. Then the energy decreases 
slowly (note that the dissipation is proportional to the surface of the sample while E is 
proportional to the volume). This lowering of energy will continue until E reaches just 
below Ec when the system falls into the absorbing state and thus the dissipation stops. 

Let us now return to the fixed energy sandpile but this time let us start from an absorbing 
state {E < Ec) and switch on the "external supervisor" who is injecting energy into the 
system. The supervisor is required to stop the injection if adding of the last energy packet 
started activity in the system. This process increases the energy E infinitesimally slowly 
and brings the system near and perhaps slightly past the threshold of activity E = Ec- 

Now, if both the dissipation at the boundaries and the "external supervisor" are present 
then the fixed energy sandpile model is nothing else but the sandpile model generating 
SOC. And we see that SOC emerges because the combined action of the dissipation and the 
"supervisor" brings the system to the critical point of the absorbing state transition of the 
fixed energy sandpile model. 

The mechanism unmasked above is rather general and present in many models of self 
organized criticality The value of recognizing this mechanism lies in making it pos- 



sible to describe and calculate scaling properties of SOC by studying "usual" nonequilibrium 
phase transitions. In particular, one may hope that field-theoretic description of SOC may 
be obtained through studies of the appropriate absorbing state transitions. 
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Of course, absorbing state transitions are numerous and it is not obvious which one is 
in the same universahty class as a given system displaying SOC. In general, continuous 
phase transitions to an absorbing state are in the universality class of directed percolation 
that can be described by the following reaction diffusion process 



A^A + A , A^O , AO^OA. (25) 

Directed percolation is rather robust to various changes in its rules but the presence of extra 
symmetries (conservation laws) may change the universality class of an absorbing state 



transition. A well known example is the parity conserving process |l54| , pSj pq] which has the 
following react ion- diffusion representation 

A ^ A + A + A , A + A^O , AO^OA. (26) 

Both of the above processes have been much investigated and the scaling properties have been 
accurately determined. Furthermore, understanding (if not complete solution) has emerged 



even on field theoretic level |^T], ^ Unfortunately, neither of the above processes have 
been directly related to models of SOC. Accordingly, the present day research is concentrated 
on absorbing state transitions which have more contact with SOC. An example is the critical 



point observed in the so called pair contact process [O] where particles diffuse only through 



the birth-death processes given by the reaction scheme 

A + A^A + A + A , A + A-^0, (27) 

where the first and the second reactions take place with probabilities p and 1—p, respectively. 
A related problem is the epidemic model ES| where the reaction scheme 



A + B B + B , B ^ A , BO^OB, (28) 

describes static healthy subjects (A) getting infected by diffusing infectious agents (B) who, 
in turn, recover with time. 

The last two models are close to the fixed energy sandpile models (and thus to SOC) 
in that both of them have an infinite number of absorbing states and their coarse-grained 
description involves an order parameter (the active particles) coupled to a static field (the 
temporarily immobile particles). It has recently been suggested that the similarity may go 



deeper, i.e. they all belong to the same universality class ||59[. This conclusion is based on 
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a field-theoretic calculation near dimension d = 6 |^ using Langevin equations which were 
suggested on phenomenological grounds for the processes p7| ) and ( pSj) At this 

point there is still a debate about both the applicability of the Langevin equations and the 
validity of the results in lower dimensions. Nevertheless, it appears that the approach of 
SOC through absorbing state transitions may be coming to an interesting and satisfactory 
conclusion. 

Of course, one should not forget that apart from the connection to SOC, absorbing state 
transitions in general constitute an important problem in the theory of NESS. The field is 
developing fast and there are many interesting details scattered across the papers. A guide 
to the models and to the extensive literature about them can be found in recent reviews 
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IV. DISTRIBUTION FUNCTIONS IN NONEQUILIBRIUM STEADY STATES 

The simplicity of the description of equilibrium system lies in the existence of the Gibbs 
distribution i.e. in the elimination of the dynamics from the calculation of averages. Al- 
though dynamics is clearly important in nonequilibrium steady states, it is not inconceivable 
that a prescription exist for a nonequilibrium equivalent of the Gibbs distribution which 
would include the essential features of the dynamics. Such a distribution function may have 
singularities as shown in simple examples [^, ^] or it may have problems with the additivity 
of the associated entropies (which is not unexpected in systems with long-range correlations) 
|6^ , |6^ . Nevertheless, a prescription with well defined restrictions on its applicability would 



be valuable and the search for nonequilibrium distribution function(s) has been going on for 
some time. 

A phenomenological approach to the above problem is the non-extensive statistical me- 
chanics [0, an approach that takes its name from the nonextensive character of the pos- 
tulated entropy. This approach has been much developed during the last decade, and not 
surprisingly, is has its success in connection with systems which have long-range interactions 



or display (multi)fractal behavior ||64 



Below we shall present a alternative approach that is somewhat less general but it is 
based on the extention of our knowledge of universality of distribution functions in strongly 
fluctuating systems. 
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A. Power laws and universality of nonequilibrium distributions 



Distribution functions of additive quantities such as e.g. the total magnetization in the 
Ising model are Gaussian in usual equilibrium systems. This Gaussianity follows from the 
central limit theorem that is applicable due to the correlations being short ranged away 
from critical points. At critical points, however, the power law correlations result in non- 
gaussian distributions. The emerging distributions are quite restricted in their possible 
shapes, however, the reason being that the distribution functions at critical points are scaling 
functions and their shape is determined by the universality class associated with the given 
critical point. 

The above observation can be used to develop a classification of nonequilibrium distri- 
bution functions. Namely, one knows that "effective" criticality (i.e. strong fluctuations 
and power-law correlations) is the norm for nonequilibrium steady states. Of course, the 
"effective" critical behavior is determined not only by the interactions but by the dynamics 
as well. Accordingly, one may expect that the scaling functions (and thus the distribution 
of macroscopic quantities) are determined by the nonequilibrium universality classes. Once 
we build a gallery of such scaling functions, we can use them in the same way as in the 
equilibrium case: We can identify symmetries and underlying mechanisms in experimental 
systems; we may find seemingly different systems belonging to the same universality class, 
and thus we can discover common underlying processes present in those systems. We can 
also use these distribution functions to find the critical dimension of a model and the ap- 
plications are restricted only by imagination. Below we shall show how to calculate these 
distribution functions in simple systems and present a few applications. 

B. Picture gallery of scaling functions 

The simplest nonequilibrium systems displaying "effective" criticality are the growing 
surfaces p7| , p8| . They are rough quite generally which means that the mean-square fluctu- 
ations of the surface diverge with system size. The roughness is defined by 



where is the area of the substrate of characteristic linear dimension L, h = ^ffh{f^ t)/^L 
is the average height of the surface, and x is a critical exponent characterizing the given 




(29) 



r 
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universality class. We shall be interested in the steady-state distribution of P{w2)dw2 and 
expect that due to criticality, the diverging scale (1^2) will be the only relevant scale in P{w2) 
and, consequently, it can be written in a scaling from 

^^(-) - (4* ((5:) ■ 

where is a scaling function characteristic of the universality class the growth process 
belong to. Below, we show how to calculate $ for a simple growth process (Edwards- 
Wilkinson equation [^) and will demonstrate that the $-s are different for growth processes 
distinct in the sense of distinct universality classes. 
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FIG. 6: Surface growth. The height of the surface above the substrate is given by h{x,t) and the 
width of the surface W2 is characterized by the mean-square fluctuation. The vertical velocity of 
the surface in general a function of the local properties of the surface, v{dxh,d'^h, ...). 

Let us begin by discussing the equations for growing surfaces. In general, deposition of 
particles on a substrate, under the assumption that the surface formed is a single valued 
function h{x,t), can be described by the equation dfh = v{h) where v{h) gives the velocity 
of advance of the interface perpendicular to the substrate (Fig.^). The velocity is usually 
written as v{h) = f -|- F{h) + t] where Vq is the average velocity due to the average rate 
of deposition, F{h) is related both to the motion of the particles on the surface and to the 
dependence of the growth on the inclination of the surface. Finally the fluctuations in the 
above processes are collected in 7] which is assumed to be a Gaussian white noise in both 
space and time. 
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A simple form for F{h) follows from the assumption that particles like to stick at points 
with large number of neighbors i.e. at large d^. Then, F{h) is approximated as F{h) = z/A^ 
and, in the frame moving with vq, one has the Edwards- Wilkinson (EW) model |]2^ of 
surface growth 

dth{x,t) = iyAh{x,t) +7]{x,t) . (31) 

This equation can be solved and one finds that the steady-state probability distribution is 
given by 

V[h{x)]= Ae-^^o(^hfdx (32) 

where a is related to u and to the amplitude of the white noise. 
Once V[h{x)] is known, P{w2) is formally obtained from 

P{w2) = {6{w2-[h^-t])) (33) 

where the average () is over all h{x) with the distribution function V[h{x)] (note that is 
a spatial average and it is still a fluctuating quantity). In practice, it is more convenient to 
calculate the generating function 

fOO _2 

G{s) = / e-'''^Piw2)dw2 = {e-'^^'-'' )) (34) 



with the above expression demonstrating why P{w2) can be calculated analytically in simple 
models. Namely, if the partition function with V can be found then the generating function 
(p^) is just the partition function of the model with a quadratic term added and such term 
usually does not spoil the solvability of the problem. Indeed, e.g. in case of the d=l EW 
model with periodic boundary conditions, the problem is reduced to the evaluation of the 
partition function of a d=l quantum oscillator thus obtaining [|6^ 
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n=l ^ / 

Now one can find the average width diverging (^2) = —dG{s) / ds\s=Q = L/{6a) in the 
L —>■ 00 limit. Using {W2) to eliminate L from (|35D , one observes that G{s) is a function 
of the product s{w2) only and, consequently the inverse Laplace transform yields P{w2) 
the scaling form (|30D . The calculation of the scaling function $(x) consist of collecting 
contributions from the poles in G{s) and one obtains 

$(x) = y ^(-1)"-^ exp ( j . (36) 

n=l ^ ' 
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Fig.|^ shows the above function displaying a characteristic shape of exponential decay 
g-TT^z/e g^|. lo^j-gQ ^ a^]2d essential singularity $(a;) ~ ^-5/2^-3/{2x) 0. 
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FIG. 7: Comparison of the scaling functions for the EW and MH models. 



On Fig.|^ we have also included the results for the so called curvature driven growth 
process which is also called the MuUins-Herring model of surface growth 1^^. This is a 
model where the rearrangement of deposited particles goes on by surface diffusion and the 
particle current jh is towards places where there are many neighboring particles i.e. Ah is 
large. This means that jh ~ VA/i and F{h) = —(A'^h. The resulting equation is called the 
MuUins-Herring (MH) equation 



dth{x, t) = -CA'^h{x, t) + ?7(x, t) . 



(37) 



The surfaces described by the MH equation belong to a universality class distinct from 
that of the EW growth. Indeed, the MH equation can be solved easily and one finds that 
{w2)mh ~ L"^ in contrast to the EW result {102) ew ~ L. Accordingly, the scaling function 
should also be different. A calculation similar to that described above for the EW case 
verifies this expectation |67] as can be observed on Fig.|^. 

An important point to remark about the comparisons of the EW and MH curves is that 
they are well distinguishable. Their maximum, their small x cutoff, and their decay at large 
X are all sufficiently different so that no ambiguity would arise when analyzing experimental 
data. Indeed, the d = 2 versions of the above $-s as well as a number of others characterizing 
various growth processes have been obtained in [B^ and it did not appear to be difficult to 



29 



pick the scaling function which was corresponding to a given set of experiments ||68[ . 

Finding out the universahty class of a growth process is one possible application if one 
has a sufficiently developed gallery of scaling functions. Below we discuss other possibilities 
for application. 

C. Upper critical dimension of the KPZ equation 



The KPZ equation is the simplest nonlinear model describing growth in terms of 
a moving interface. It differs from the EW model by taking into account that the surface 
grows in the direction of its normal provided the incoming particles have no anisotropy in 
their arrival direction. Then the z component of the velocity of the surface has a correction 
term proportional to (V/i)^ as shown on Fig.^ and the equation in lowest order in the 
nonlinear ities becomes the so called KPZ equation 

9t/i = i/V^/i + A(V/i)^ + r7. (38) 

Here v and A are parameters, while riif^ t) is again a Gaussian white noise. The steady state 
surfaces generated by (^) appear to be rough (critical) in any dimension if the nonequi- 
librium drive (A) is large enough. Since eq. (|38D gives account of a number of interesting 
phenomena (Burgers turbulence, directed polymers in random media, etc.) lots of efforts 



have been spent on finding and understanding the scaling properties of its solutions |28 
Nevertheless, a number of unsolved issues remain, the question of upper critical dimension 
(du) being the most controversial one. On one hand mode-coupling and other phenomeno- 



logical theories suggest that (i^ = 4 |^ while all the numerical work fail to find a finite du 



and the indication is that du = oo ||7^ . Below I would like to show how the scaling functions 



of the roughness can shed some light on this controversy . 

Let us begin with the observation that scaling functions do not change above du- Thus 
if we build $(x) in dimensions d = 1 — 5 and observe that they differ significantly in d = 4 
and 5 then we can conclude that du > 4. Since $(a;) cannot be exactly calculated for d >2 
we must evaluate it through simulations with the results displayed on Fig.^]. 

As one can see on Fig.||, the scaling functions change smoothly with d. The $(x)-s 
get narrower and more centered on x = 1 with increasing d, and there is no break in this 
behavior at = 4. The equality of the d=4 and 5 scaling functions appears to be excluded. 
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FIG. 8: Roughness distribution for KPZ steady-state surfaces in dimensions d = 1 — 5. 

Of course, our result is coming from numerical work with the same general conclusion 
as in previous studies. So, why is it more believable? Because one of the main criticism of 
numerical studies does not apply to it. Namely, no fitting parameters and fitting procedures 
are used in contrast to the usual determination of critical exponents. One just builds the 
histograms, calculates the averages to determine the scaling variable and plots the scaled 
histogram. This is clear and well understood but there is another remarkable feature in 
these scaling functions the origin of which is less obvious. Namely, the scaling functions in 
principle should depend on the size of the system. 



(39) 



What is observed, however, is that the L dependence is practically all in (^2)^ and L 
dependence of the shape of $ (explicit dependence on L) disappears already at small L. 
These are important points and the KPZ application of the scaling functions was mainly 
chosen to emphasize them. 

It seems that these scaling functions are versatile tools which can be used in computer 



science |^ as well as in understanding the propagation of chemical fronts [^. An interesting 
application was e.g. the establishment of a connection between the energy fluctuations in a 
turbulence experiment and the interface fluctuations in the d = 2 Edwards- Wilkinson model 



74 1 , and thus prompting a search for an interface interpretation of the dissipative structure 
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in the turbulent system [^. In another case, it helped to make a link between the much 
studied 1// noise and the extreme value statistics [^. Since the effective criticality is a 
real feature of many nonequilibrium system, we expect that many more use will be found 
for the scaling functions discussed in this section. 

V. QUANTUM PHASE TRANSITIONS 

Quantum critical points are associated with the change of the symmetry of the ground 
state of a quantum system as the interactions or an external field (control parameter) are 
varied. They have been much investigated in recent years with the motivation coming from 
solid state physics [^]. Namely, the strongly-correlated electron systems often produce 



power law correlations, and the origin of the observed scale- invar iance is suggested to be 
the presence of a quantum critical point at T = provided the effect of the quantum phase 
transition is felt at finite T as well. 

From our point of view, the quantum phase transitions are interesting because they 
are good candidates for studying the effects of a nonequilibrium drive on well established 
symmetry-breaking transitions. The advantage of these systems is that there is no arbitrari- 
ness in their dynamics (it is given by quantum mechanics), the one-dimensional systems are 
simple with examples of exactly solvable models displaying genuine critical phase transi- 
tions (see e.g. the transverse Ising chains discussed below) and, furthermore, there is much 
previous work to build on. 

The only problem is how to force a quantum system into a non-equilibrium steady state. 
An obvious way is to attach two heat baths of different temperatures at the two ends of a 
spin chain. Unfortunately, this makes the problem unsolvable (even numerically) for any 
reasonable size system and thus it is practically impossible e.g. to draw conclusions 
about the long-range correlations generated in the system. Below we show a way to avoid 
the problem of heat baths. The idea is that the nonequilibrium steady states always carry 
some flux (of energy, particle, momentum, etc.). Thus a steady state that is presumably not 
very far from the one generated by boundary conditions may be constructed by constraining 
the quantum system to have a flux equal to the one generated by the boundary conditions. 
For example, in the case of the transverse Ising chain treated below, we shall constrain the 
system to carry an energy current and will investigate the correlations in this constrained 
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state. ^ 



A. Spin chains with fluxes 

As a simple model with critical phase transition, we consider the d = 1 Ising model in a 
transverse field h which has the following Hamiltonian: 

Hi = -Yl [''i^i^i + ifl) ■ (40) 

Here the spins cr" [a = x,y,z) are represented by 1/2 times the Pauli matrices located at 
the sites i = 1,2,..., N of a one-dimensional periodic chain (c^+i = cr"). The transverse 
field, h, is measured in units of the Ising coupling, J, which is set to J = 1 in the following. 

This model can be solved exactly |]8^ , ^ and it is known that a second order phase 
transition takes place in the system as h is decreased. The order parameter is the expectation 
value (ax) i.e. (cr^) = for /i > 1, while (ax) 7^ for /i < 1 and he = 1 is a critical point. 
The scaling behavior at and near belongs to the d = 2 Ising universality class. 

In order to constrain the above system to carry a given energy fiux Je we shall use the 
Lagrange multiplier method. Namely, we add a term XJe to the Hamiltonian where Je is 
the local energy fiux operator summed over all sites, and find that value of A which produces 
a ground state with the expectation value (Je) = Je- 

The above scheme requires the knowledge of the local energy current, Je. It can be 
obtained using the quantum mechanical equation of motion for the energy density ii = 
i/h[Hj, Ee], and representing the result as a divergence of the energy current ii = Je — Je+i- 
The calculation yields (/i = 1 is used in the following) 

Je = \<yJ{<yU - ^hi) (41) 

and this allows to construct the 'macroscopic' current Je = Yle ^i-- Adding it to Hi with a 
Lagrange multiplier, —A, 

H = Hi-XJe . (42) 

we obtain the Hamiltonian whose ground states with (Je) = 7^ will give us information 
about the current carrying states of Hj. 

^ It should be noted that this section was not discussed during the main lectures of the course. It was 
described only in a seminar for interested students. 
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In order to avoid confusion, we emphasize that the energy current, Je-, is associated with 
Hi and not with the new Hamiltonian, H. We also note that H is just another equihbrium 
Hamiltonian, it differs from Hi by an extra term which breaks the left-right symmetry of 
Hj. Finding the ground state of if, however, gives us the minimum energy state of Hi which 
carries an energy current, Je = (Je)- Thus the ground-state properties of H provide us 
with the properties of the nonequilibrium steady states of the transverse Ising model. 

It turns out that [Hj, Je] = and H can be diagonalized by the same transformations 
which diagonalize Hi ||8^, and one arrives to a system of free fermions with a spectrum of 
excitation energies given hj Ug = \Aq\ where 

= i^l + h'^ + 2hcosq + ^sing . (43) 

with the wave numbers restricted to — vr < g < vr in the thermodynamic limit (A^ —>■ oo). 

1.5 
1 

< 0.5 


-0.5 

-71 7C 

q 

FIG. 9: Spectrum of the transverse Ising model in the presence of a field (A) which drives the 
current of energy. The excitation energies are given as ujg = \Aq\. Increasing the drive makes the 
ground-state change at a critical A = Ac (Ac = 3 for h = 2/3) when negative energy states start to 
appear and get occupied. The qualitative picture is the same at all transverse fields h. 

Fig.p displays the spectrum for h = 2/3 and various A and one can see that the g —>■ —q 
symmetry of the spectrum is broken for A 7^ 0. Nevertheless, for small A the ground-state 
remains that of the transverse Ising model (A = 0) since Ag > and the occupation number 
representation of the ground state does not change. Accordingly, no energy current flows 
( Je = 0) for A < Ac. This rigidity of the ground state against the symmetry-breaking field 
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q>-~-__.,-'-'''q+ 





which drives the energy current is a consequence of the facts that the fermionic spectrum of 
the transverse Ising model has a gap and that the operator Je commutes with Hj (similar 
rigidity is observed in the studies of energy flux through transverse XX chain [B^ ) . 
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FIG. 10: Phase diagram of the driven transverse Ising model in the h — X plane where h is the 
transverse field while A is the effective field which drives the flux of energy. Power-law correlations 
are present in the nonequilibrium phase (Je ^ 0) and on the Ising critical line in the equilibrium 
phase {Je = 0, dashed line). 

The ground-state properties do change when < in an interval [g_,g_|.] and these q 
states become occupied. Due to the resulting asymmetry in the occupation of the q and —q 
states, the energy current becomes nonzero. The line Xc{h) which separates the region of 
unchanged transverse Ising behavior from the Je ^ ^ region is obtained from the conditions 



Aq = and dAg/dq = 0, and is displayed on the phase diagram (Fig.|T^ as a solid line. 
Another phase boundary on Fig.|lO| is shown by dashed line. It separates the magnetically 
ordered {h < 1, X < 2/h) and disordered {h > 1, X < 2) transverse Ising regions. Since the 
ground state is independent of A for A < Ac, one has the same second order transition across 
the dashed line as at h = 1 and A = i.e. it belongs to the d = 2 Ising universality class 



Clearly, one can view the region A < Ac as an equilibrium phase while the A > Ac region 
as a nonequilibrium one since there is a nonzero energy flux through the latter. This flux 
can actually be calculated easily with the simple result 



Je 



{Je/N) = (47r)-V(l-4/A^) {h'-^/X') 



(44) 
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Apart from the fact that 7^ 0, the A > Ac region should also be considered as a 
distinct phase since the long-range magnetic order existing for h < 1 breaks down when 
Je ^ and the magnetic correlations become oscillatory with amplitudes decaying as a 
power of distance. Indeed, this can be seen by investigating the (erf cr|i(_„) correlations which 
can be expressed through Pfaffians [BT| and thus making possible numerical calculations for 
n < 100. In the presence of long-range order one should have {<J^<7^^,J (<^£ )^ 7^ for 
n ^ 00 while we find that the correlations decay to zero at large distances as 

K^.+n)~^^cos(A;n) (45) 



where the wavenumber, k = arccos (2/A/i). The above result (|45|) is coming from numerics 
and it is exact in the A ^ 00 limit where the correlations are those of the d = 1 XX model 



m 



One can observe power-law correlations for A > Ac in other physical quantities as well. 
For example, the envelopes of both (o"|cr|^„) and {JeJe+n) correlations behave as in 



the large n limit I^Oj. Thus we arrive to the main conclusion of this section, namely that a 
simple, exactly soluble quantum system shows power-law correlations in the current carrying 
state in agreement with the notion that power-law correlations are a ubiquitous feature of 
nonequilibrium steady states. 

Actually, remembering that power-law correlations in quantum models are associated 
with a gapless excitation spectrum, we can reformulate the transverse Ising model result to 
see a general connection between the emergence of power-law correlations and the presence 
of a current. Indeed, let us assume that a system with Hamiltonian Hq has a spectrum 
with a gap between the ground-state and the lowest excited state. Furthermore, let J be 
a 'macroscopic' current of a conserved quantity such that [Hq, J] = 0. Generally, there is 
no current in the ground state and adding — AJ to Hq does not change the (J) = result 
for small A. Current can flow only if some excited states mix with the ground state and, 
consequently, a branch of the excitation spectrum must come down and intersect the ground- 
state energy in order to have (J) 7^ 0. Once this happens, however, the gap disappears and 
one can expect power-law correlations in the current- carrying state. Admittedly, the above 
argument is not strict and is just a reformulation (in general terms) of what we learned 
from the transverse Ising model. We believe, however, that the above picture is robust and 
suggestive enough to try to find other soluble examples displaying the fiux power-law- 
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correlations relationship. 



B. Quantum effective interactions 

When using the Lagrange multiplier method, one assumes that the flux generated by 
boundary conditions can be replaced by the effective interactions contained in an appropri- 
ately chosen global flux J. These interactions are generally short ranged since the flux is 
usually a sum of local terms. The short-range nature of the effective interactions is actu- 
ally not in contradiction with the power-law correlation being generated since the Lagrange 
multiplier gets tuned in order to achieve a given flux and, furthermore, the tuning is not 
quite trivial since A must be increased past a critical value in order to have e.g. a nonzero 
flux energy. The question nevertheless arises whether adding a flux term in H was an ade- 
quate description for the nonequilibrium steady state which is expected to display power-law 
correlations. 

In order to investigate the above question one would have to solve the problem with the 
boundary drive but, as discussed above, an exact solution does not seem to be feasible. 
Instead, however, one can prepare initial conditions which will lead in the long-time limit to 
a steady flux. Then the steady state obtained in this natural way can be compared to the 
one found by the Lagrange multiplier method. This program has been carried out |^ for 
the XX chain defined by the following Hamiltonian 

N 

Hxx = -Y. i^i^i+i + + h^i) ■ (46) 

1=1 

In this model, the transverse magnetization = Y2i '^t conserved and one can investigate 
the nonequilibrium states which carry a given magnetization flux by using the Lagrange 
multiplier method At the same time, the model is simple enough so that one can solve 
the time dependent problem where a steady magnetization flux is achieved by starting with 
an inhomogeneous initial state that is the ground state at fixed magnetization but with 
m = (s^) reversed from rriQ for n < to —tuq for n > 0. The time-evolution of this step-like 
initial state can be followed exactly and the magnetization profile emerging in the long-time 



limit is shown on Fig.^. The remarkable feature of this magnetization profile is the middle 
part which is an m = homogeneous state carrying a magnetization flux j(mo). Comparing 
this state to the one generated by adding the flux term to the Hamiltonian and fixing the 
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FIG. 11: Time evolution of the magnetization profile starting from a step-like initial conditions 
shown as solid line. There are two fronts going out to ibco. They diminish the magnetization 
and leave behind a homogeneous (s^) = state. In the scaling limit t oo, n/t — > x, the 
magnetization m{n,t) ~ ^{n/t) is given by <^{x) = mo for —1 < x, <I>(x) = mo — vr^^ arccos(a;) for 
— 1 < X < — cos(-7rmo), ^{x) = for — cos(7rmo) < a; < 0, and ^{x) = — <!•(— x) for a; > (dashed 
line). 

Lagrange multiplier to have the same j{mo), we find that various expectation values such 
as the energy, the occupation number in fermionic representation are all equal in the two 
states. Thus the Lagrange multiplier yields a correct description of the states carrying a 
magnetization flux. 

It should be noted, however, that recent calculation with an inhomogeneous initial state 
of different temperatures (Ti for a; < and T2 for x > 0) has yielded a different result for 
the asymptotic state carrying an energy finx |^6[|. Namely, it was shown that, at least in the 
neighborhood of a; = and in the t — > 00 limit, the properties of the flux-carrying state can 
be interpreted in terms of the ground state of an effective Hamiltonian 

N N 
n=l j=l 

where Qj"'* is a product of local operators at sites j and j + n, and the interaction is of 
long-range type since /x„ ~ 1/n |^6| (remarkably, the first two operators Q^p and Qj^^ are 
those appearing in the Lagrange multiplier treatment of the energy flux in the XX chain). 
Although the homogeneity of the asymptotic state was not shown and thus the comparison 
may be questioned, the above result indicates that the Lagrange multiplier approach may 
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-m, 







be only a first approximation in describing the fiux-carrying states. 

In summary, the studies of quantum systems described above strengthen the view that 
fluxes generate long-range correlations. Furthermore, the quantum systems also give a simple 
picture of how the emergence of these correlations is related to the closing a gap in the 
excitation spectrum above the ground state. 

VI. OUTLOOK 

There are topics which are important but were not discussed in these lectures. To mention 
a few, there is a large body of work on one-dimensional systems displaying nonequilibrium 
phase transitions, on orderings of granular gases under shear, on pattern formation with 
phase transitions described by the complex-coefficient Landau- Ginzburg equation, and the 
list could be continued. My choice of topics mainly reflects my past work and my attempts 
to develop simple starting points for making inroads into the beautiful but rather difficult 
field of far from equilibrium phenomena. 

Finally, I was asked to provide entertainment for readers by trying to guess the future 
developments in connection with nonequilibrium orderings. Well, one of the present problem 
of the field is the lack of simple experimental systems which can be be taken far enough 
from equilibrium and compared to elementary models of NESS. Search for such systems 
will intensify and I expect that there will be a shift towards biological problems. There the 
condition of being far from equilibrium is satisfied and there may be surprisingly simple 
phenomena under the guise of complicated pictures. This line of research may in the future 
meet up with game theories generalized to take into account spatial structures. 

Search for better understanding of the emerging effective interactions will also continue, 
just as the sorting out of the absorbing-state transitions (surprising connections may be still 
found there, in addition the existing one to SOC). I also believe that theory of nonequilibrium 
distributions will be much developed, and limiting distributions such as the ones emerging 
in extreme statistics will have a much wider use in physics. These are more or less safe bets. 
And then there is the unpredictable part of future. 
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